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Direct numerical simulation of homogeneous nucleation and growth in a phase-field 

model using cell dynamics method 

Masao IwamatsiJ^ 
Department of Physics, General Education Center, 
Musashi Institute of Technology, Setagaya-ku, Tokyo 158-8557, Japan 
(Dated: February 8, 2008) 

Homogeneous nucleation and growth in a simplest two-dimensional phase field model is numeri- 
cally studied using the cell dynamics method. Whole process from nucleation to growth is simulated 
and is shown to follow closely the Kolmogorov-Johnson-Mehl-Avrami (KJMA) scenario of phase 
transformation. Specifically the time evolution of the volume fraction of new stable phase is found 
to follow closely the KJMA formula. By fitting the KJMA formula directly to the simulation data, 
not only the Avrami exponent but the magnitude of nucleation rate and, in particular, of incubation 
time are quantitatively studied. The modified Avrami plot is also used to verify the derived KJMA 
parameters. It is found that the Avrami exponent is close to the ideal theoretical value m = 3. The 
temperature dependence of nucleation rate follows the activation-type behavior expected from the 
classical nucleation theory. On the other hand, the temperature dependence of incubation time does 
not follow the exponential activation-type behavior. Rather the incubation time is inversely pro- 
portional to the temperature predicted from the theory of Shneidman and Weinberg [J. Non-Cryst. 
Solids 160, 89 (1993)]. A need to restrict thermal noise in simulation to deduce correct Avrami 
exponent is also discussed. 

PACS numbers: 64.60.Qb, 68.18.Jk, 81.10.Aj 
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I. INTRODUCTION 



The dynamics of phase transformation by nucleation 
and growth of a new stable phase from a metastable 
phase is a very old problern, which has been studied for 
more than half a century [J, [3, Q from a fundamental 
point of view as well as from technological interests. Ex- 
perimentally 0, Q, the dynamics of phase transforma- 
tion has been believed to follow the classical Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) picture of nucleation and 
growth ^,[1,0, [I]- According to their picture, the phase 
transformation proceeds via the nucleation and subse- 
quent growth of nuclei by the interface-limited growth. 
Therefore, the classical nucleation theory and the steady- 
state growth are assumed. Their picture is integrated 
into the well-known KJMA formula for the fraction of 
transformed volume. 

Theoretically, on the other hand, the various variant 
of phase-field [3, [13, [HI model has been routinely used 
to study the dynamics of phase transformation. Essen- 
tially similar models called Cahn-Hilliard mode l |13 | and 
time-dependent Ginzburg-Landau model [l^ Il4| have 
also been used extensively. It is not obvious, however, 
that the dynamics of phase-field model is in accord with 
the KJMA picture. In particular, the time scale of phase 
transformation including the transient nucleation has not 
been fully studied. For more than a decade ago. Vails 
and Mazenko lj| studied the dynamics of phase trans- 
formation in a phase-field model. They studied nucle- 
ation and growth rather phenomenological way and did 



not pay attention to the connection to KJMA formula. 
Very recently, Jou and Lusk 9], and Castro fld\ stud- 
ied the connection between phase-field model and KJMA 
formula. However, the former introduced critical nucleus 
(seed of new phase) artificially and the latter paid more 
attention to the heterogeneous nucleation. Therefore, the 
nucleation process is artificially controlled in their stud- 
ies d, [13] ■ Very recently, Granasy et al. |15| extensively 
studied the validity of the KJMA picture in the phase- 
field model of spherulite. They paid more attention to 
the Avrami exponent rather than the time scale of phase 
transition. Therefore, virtually there has been no de- 
tailed study of the time scale of phase transition including 
both the nucleation and the growth in phase-field model 
on the same footing. 

In this paper, we will use the phase-field model with 
thermal noise to study the homogeneous nucleation and 
the growth in a unified manner. In Sec. II we present a 
short review of the classical KJMA picture of nucleation 
and growth. In Sec. Ill, we present the phase-field model 
and the cell dynamics method. Section IV is devoted to 
the results of numerical simulations. Finally Sec. V is 
devoted to the conclusion. 



II. CLASSICAL KJMA PICTURE OF 
NUCLEATION AND GROWTH 

The time evolution of the volume fraction / of trans- 
formed volume for two-dimensional system is predicted 
from KJMA (Kolmogorov-Johnson-Mehl-Avrami) [l|, 0, 
0, H, [ill theory as 
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/ = l-exp [-^n2{t)], 



(1) 



2 



with 



n2{t) = / I{t')R\t;t')dt\ 
Jo 



(2) 



where I{t) is the (time-dependent) nucleation rate and 
R{t; t') is the radius of nucleus at time t that was nucle- 
ated at t'. 

These equations are further simplified if [l6j 

(i) the interfacial velocity Vs is time- and size- 
independent, and 

(ii) the radius of critical nucleus is infinitesimally 
small, and 

(iii) the nucleation rate (for critical nucleus) I{t) is ap- 
proximated by the (time-independent) steady state 
nucleation rate Is, 

then we have a linear time dependence of radius of nuclei 
R{t; t') = Vs{t — t'), and the integral in Eq. ^ can be cal- 
culated analytically to give the classical KJMA formula: 



/ = 1 - cxp 



(3) 



However, if we take into account the fact that the 
radius of critical nucleus i?* is finite, we have to shift 
the origin of time scale toward the past by the amount 
to — R^,/vs. Then, we have to replace time t by [iTj 



(4) 



According to the theory of Shneidman and Weinberg [T^| , 
on the other hand, if we take into account the transient 
nucleation rate I{t) as well as the size dependence of 
interfacial velocity v, we have to replace the origin of 
time t toward future by 



t t 



1 



to 
t 



In 



t 

to k^T' 



(5) 



where W* is the energy barrier to form the critical nu- 
cleus and T is the absolute temperature, since we have 
to wait for the nucleation rate as well as the interfacial 
velocity to reach the steady state values. Therefore two 
contributions Eq. ^ and Eq. ^ have opposite signs and 
we have to replace time t in Eq. ([3]) by 



t ^ t-ti, 



(6) 



where iinc is usually called " incubation time" , which now 
consists of two contributions Eq. ^ and Eq. 



tu 



to In 



t 

to ksT 



to. 



(7) 



where the time t ^ to is the timescale of the total process 
of phase transformation. At low temperature (small T) 
or low undercooling (small e and, therefore, large W,), 
critical nucleus is rarely formed and the time necessary 



to reach the steady state would be long. Then the first 
term of Eq. ([7]) dominates and we would have a posi- 
tive tine > 0. Usual condition of the classical nucleation 
theory W^/ksT ~ 40 - 70 satisfies this condition. 

Furthermore, if we have anomalous power-law depen- 
dence of nucleation rate I{t) ~ /pt" and interfacial ve- 
locity v{t) ~ vot^ , then we have a generalized KJMA 
formula 



/ = 1 - exp 



t tin 



where 



m = 3 + a + 2/3 



(8) 



(9) 



is the so-called Avrami exponent, ti^c the incubation 
time, and tg^ is the growth time, which is determined 
from the nucleation rate /q and the interfacial velocity 
Vo through 



1 



i-gr 



— ^0«0 
TO 



(10) 



The ideal case Eq. (jS]) is given by the formula with a — 
P = and TO = 3. 

In the classical nucleation theory (CNT) [Ell, [1], the 
steady-state nucleation rate ~ /q, that is the number 
of critical nuclei which appear per unit time and unit 
volume, is usually given by the activation form 



Is oc exp — 



k^ 



(11) 



where is the nucleation barrier of critical nucleus that 
appeared in Eq. ((5]). Then, the growth time is given by 
the activation form: 



tgr cx exp 



\mk^T 



(12) 



Similarly, by replacing the time scale t in Eq. ([7]) by t^ 
in Eq. (fT2|) . we have 



to / mkBT 



(13) 



The two time-scales of phase transformation, the incu- 
bation time tine and the growth time tgr will increase 
as the temperature T is lowered. But the temperature 
dependence is different. The growth time increases expo- 
nentially and it increases much faster than the incubation 
time as the temperature is lowered. 

The classical expression of the energy barrier for 
the two-dimensional circular nucleus is given by 



'Ah 



with radius of critical nucleus R^: given by 



(14) 



(15) 



3 



where a is the interfacial energy of nucleus and Aft is the 
free energy difference between the metastable and the 
stable phase. 

In section IV we will check if this classical KJMA pic- 
ture of nucleation and growth PjSHjBIHIj in particular 
Eqs. (HI), (fT2|) and (fTS]) are vahd in the phase- field model 
using cell dynamics method. A similar study to test 
the validity of KJMA picture in Ising-type spin model 
was conducted by Shneideman et al. [l7[ and Ramos et 
al. [11. 



III. PHASE-FIELD MODEL AND CELL 
DYNAMICS METHOD 

In the phase-field model, the dynamics of phase trans- 
formation is described by the standard isothermal phase- 
field equation [13, [HI, [S] for continuous variables 



dtp 
'dt 



SF_ 



(16) 



where S denotes the functional differentiation, -0 is the 
non-conserved order parameter, and F is the free energy 
functional. This free energy is usually written as the 
square-gradient form: 



[D{\/tP)^ + h{ij)] dr. 



(17) 



The local part h{tp) of the free energy F determines the 
bulk phase diagram and the value of the order parameter 
in equilibrium phases. 

In the cell dynamics method, the partial differential 
equation (jl6p is replaced by a finite difference equation 
in space and time in the form 



ip{t + l,n) = M[iP{t,n)], 



(18) 



where the time t is discrete and an integer, and the space 
is also discrete and is expressed by the integral site index 
n. The mapping M is given by 

M[i;{t, n)] = -/(V'(i, n)) + D [« ^{t, n) » -^{t, n)] , 

(19) 

where /(V') — dh{ip)/dip, and the definition of • • • ^ 
for the two-dimensional square grid is given in [Tol . [20j . 
We use the map function /(f/') directly obtained from the 
free energy ft(V') [HI [11| instead of the standard tanh 
form originally used by Oono and Puri [l^, [2^, which 
is essential for studying the subtle nature of nucleation 
and growth when one phase is metastable and another is 
stable. 

The local part of the free energy /i(?A) we use is [ol. [2l| 



/i(V) = -V'(i - V-)' + 



T 



(20) 



This free energy is shown in Fig. [1] where one phase at 
■0 = is metastable while another phase at ■)/) = 1 is 



stable. The free energy difference Ah between the stable 
phase aX ip = 1 and the metastable phase at V' = is 
determined from the parameter e: 



Ah = ft(0 = 0) - h^iP 



(21) 



We will use the terminology undercooling to represent e. 
The metastable phase at ■0 = becomes unstable when 
e = 1/3, which defines the spinodal. 
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FIG. 1: Model double- well free energy h{'ip) defined by 
Eq. (|2Up . The undercooling parameter e determines the free 
energy difference Ah = e/4 between the depth of two wells. 
The phase with ip — is metastable while the one with tp — 1 
is stable. 

The interfacial energy a in Eqs. and can be 
estimated from the interfacial profile calculated from the 
Ginzburg-Landau equation SF/Sip = for the order pa- 
rameter tp at the two-phase coexistence (e = 0) [23j . 
Then the e-dependence of the surface tension a could 
be neglected, and the energy barrier and the radius 
of critical nucleus in Eqs. and becomes inversely 
proportional to the under cooling e: 

W^oc-, R^oc- (22) 
e e 

These formulae are incorrect when the spinodal is ap- 
proached when e 1/3 as the surface tension will van- 
ishes (cr — + 0) at the spinodal. Then, both the energy 
barrier and the radius i?* of critical nucleus vanish 
at the spinodal. Equation ([^ cannot be used near the 
spinodal. 

In the previous paper [2l| , using this phase- field model 
and the cell dynamics method, we could successfully sim- 
ulate the interface-limited growth of a single nucleus [2^ 
with a constant growth velocity v which is close to the 
theoretical prediction [2l|, 



(23) 



1 

V — — -, , 
2 V 2 

The interfacial velocity v is constant during evolution 
and does not depend on the temperature, that is in ac- 
cord with the assumption (i) of the original KJMA for- 
mula Eq. ([3]). Therefore, the temperature dependence of 



4 



growth time tg^ in Eq. ([TO)) comes solely from the temper- 
ature dependence of nucleation rate Ig given by Eg. (|TT|). 

By introducing the critical nucleus artificially [3, [2l| , 
the growth of the multiple of nuclei with the site- 
saturation condition as well as with the continuous ho- 
mogeneous nucleation condition could also be studied. 
In this case, the critical nucleus is not introduced in 
the region where the phase transformation has already 
started. Introduction of an extra nucleus on the bound- 
ary of phase transformed region will accelerate the phase 
transformation and will result in the increase of apparent 
nucleation rate. With such a care, the phase-field model 
can reproduce the theoretically predicted ideal Avrami 
exponent m ~ 2 for the site saturation case and to ~ 3 
for the continuous nucleation case '2l'|. 

In order to simulate not only the growth process but 
the nucleation process in the phase field model, we add 
thermal noise in this paper in Eq. (|18p : 

^p{t + l,n) ^ M[^j{t,n)]+^{t,n) (24) 

The thermal noise ^(t) should be added to only those re- 
gions where the phase transformation has not yet started. 
This process is simulated by adding thermal noise ^{t, n) 
only when the phase field is smaller than some cutoff 
value ipc- Therefore stochastic Eq. ([24]) is used only when 
tp{t, n) < ipc otherwise deterministic Eq. (fT8|) is used. We 
set ipc = 0.2 thought this paper. Too small V'c < 0.2 sup- 
press nucleation completely, while larger tpc > 0.2 will en- 
hance the nucleation rate and will lead to larger Avrami 
exponent m > 3 and non-linear behavior of nucleation 
and growth parameters. The limitation of this ad hock 
method will be discussed later in section [TVl 

The thermal noise £^{t, n) is related to the absolute tem- 
perature T from the fluctuation-dissipation theorem as 

{at,ri)C{t',n')) ^kBTS„,n'St,t'- (25) 

In this paper, we will use a uniform random num- 
ber rangin g fr om —\/t to +\/t for the thermal noise 
^{t,n) [l9|, |20|. Then, the parameter r is proportional 
to the absolute temperature: 

T cx T, (26) 

and the temperature T is included through the thermal 
noise r. 

IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

Using the cell dynamics method presented in the previ- 
ous section we have simulated the nucleation and growth 
in a two-dimensional phase field model. The system size 
is fixed to 512 x 512. Throughout this work, we set 
D = l/2 and ijjc = 0.2. 

Figure [5] shows a typical pattern of evolution of a new 
phase (white) in the metastable old phase (black). The 



initial phase is an unstable phase with the order pa- 
rameter i/} = (black). In contrast to the previous 
method where new nuclear embryos are artificially in- 
troduced near circular nuclei of new phase ap- 
pear spontaneously from the thermal noise r and start 
to grow (Fig. [2]). As the radius i?* of critical nucleus is in- 
versely proportional to the undercooling e from Eq. (|15p , 
the growing nuclei are larger for low undercooling e = 0.2 
(Fig. ma)) than for high undercooling e = 0.3 (Fig.[2l^b)). 
The time scale of phase transformation t-mc and igr are 
much shorter for high undercooling e = 0.3 near the spin- 
odal e = 1/3 than for low undercooling e = 0.2. 
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FIG. 2: A typical evolution pattern of the phase transforma- 
tion for 512x512 system when (a) the undercooling e = 0.2 
and (b) e — 0.3 at a low temperature 1/r = 40. Note that 
the spinodal is at e = 1/3. 

Figure [3] shows a typical shape of the time evolution 
of the volume fraction / of transformed volume for two- 
dimensional 512 X 512 system averaged over 10 samples 
with different sequences of random number when e — 0.25 
and = 25. Since the statistical standard deviations 
of 10 samples are too small to be visible on the scale 
of the figure, we did not show the error bars in figure 
m The theoretical curve (solid line) is obtained by the 
least-square fitting of three parameters, igr, iinc and m 
in Eq. ([8]) to the simulation data. Only those simulation 
data within 0.05 < / < 0.095 [11] is used for least-square 
fitting to avoid numerical error. The values for these 
three parameters differ less than 1% for larger 1024 x 1024 
system and at most 5% for smaller 128 x 128 system. 
Therefore the use of rather small 512 x 512 system is 
justified. In contrast to the previous studies 0, [2l|, it 
can be seen from Fig. [3] that the KJMA formula Eq. ([5]) 
can reproduce the growth curve of / almost perfectly. 

The Avrami exponent to, the growth time igr and the 
incubation time iinc determined for various temperatures 
T and undercooling e are summarized in Table [H The 
Avrami exponents are almost all around the ideal value 
3. Therefore, the classical KJMA picture of nucleation 
and growth can successfully describe the phase transfor- 
mation kinetics in our simple phase field model. 
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FIG. 3; The time evolution of the volume fraction / average 
over 10 samples for 512 x 512 system (crosses) compared with 
the theoretical prediction Eq. ([8]) (solid lines) when e = 0.25. 



TABLE I: The Avrami exponent m, the growth time tgr, 
and the incubation time tine determined from direct fitting 
of Eq.® to simulation data and those determined from the 
modified Avrami plot Eq. 
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Figure SJa) shows the Avrami plot corresponding to 
Fig. [3l From Eq. ([8]) we expect a linear relation 

In [- ln(l - /)] = m In (t - tine) + constant. (27) 

between In [— ln(l — /)] and In (< — ii„c) which is the so- 
called Avrami plot. The slope gives the Avrami expo- 
nent m = 2.8005 ~ 2.80 in Fig. gfa) that is very close 
to m = 2.83 for e = 0.25 and t'^ = 25 in Table U de- 
duced directly by fitting the KJMA formula Eq. dS]). In 
contrast to the previous approach where the critical nu- 



cleus is artificially introduced, [l,[2l[, our simulation data 
in Fig. m^a) almost perfectly fits to the KJMA formula 
Eq. (P7|) . This is partly due to the fact that the incu- 
bation time is unambiguously defined and determined in 
our simulation, while there is uncertainty of the incuba- 
tion time in previous works [21] as the nucleus of finite 
size is artificial introduced during the evolution. The ex- 
ponents m for other e and r deduced from the slope of 
the Avrami plot are almost the same as those obtained by 
direct fitting of Eq. ([5]) to the simulation data tabulated 
in Table [J The results differ at most 10%. 

It should be noted that the Avrami plot in Fig. [H^a) 
cannot be used directly to deduce the Avrami exponent 
m as the plot needs accurate pre-determination of the 
incubation time tine- Inappropriate choice of this in- 
cubation time is known to lead to non-linear Avrami 
plot [2l[. In order to remedy this deficiency, Mao and 
Altounian [26l| proposed modified Avrami plot. By dif- 
ferentiating Eq. ([5]), we have 



dt 



m 
t7 



which is combined with Eq. 



(1-/): 

to give 



(28) 



ln[-ln(l-/)] 



m 



771—1 



■In 



771 — 1 



■In 



df/dt 



1-/ 
_ ■ (29) 
Eq. ((29)) is called modified Avrami plot [20|. Therefore, 
we expect a linear relation between In [— ln(l — /)] and 
ln[df /dt/ {1 ~ f)]. Now, the predetermination of incu- 
bation time tine is unnecessary, and the Avrami expo- 
nent can be deduced from the slope m/(m — 1), and 
the growth time <gr can be deduced from the constant 

— (to/ (?77 — 1)) In (TO/tgr). 

This modified Avrami plot in FigurelU^b) clearly shows 
a predicted linear relation. Since we use cell dynam- 
ics method, we have replaced the differentiation df/dt 
by the finite difference in Eq. (|29p . The Avrami ex- 
ponent deduced from this modified Avrami plot to = 
1.53441/(1.53441 - 1) ~ 2.87 is very close to to = 2.83 
for e = 0.25 and = 25 in Table U deduced directly 
by fitting KJMA formula Eq. ([S]) to simulation data. 
Therefore, the validity of KJMA picture in our phase 
field model is further confirmed. The Avrami exponents 
TO and the growth times tgr deduced from this modified 
Avrami plot differs at most a few percent from those de- 
duced from direct fitting (Table [l|. 

Although, the validity of KJMA picture represented 
by the KJMA formula in Eq. ([5]) in our phase field model 
is confirmed, the Avrami exponent to obtained by fitting 
to the simulation data shows slight deviation from the 
ideal value m = 3 as shown in Table [H which means that 
the time dependence of nucleation rate / ~ lot"' and 
interfacial velocity v ~ vot^ with a Q or /3 ^ exist. 

Figure[5]shows the incubation time tine estimated from 
the least-square fitting of KJMA curve in Eq. ^ to the 
simulation data and tabulated in Table U as the function 
of temperature parameter l/r for various undercooling 
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FIG. 4: (a) The Avrami plot of / (crosses) of Figure [3] com- 
pared with theoretical predictions Eq.([8)l (solid lines). The 
incubation time tine obtained from the least-square fitting in 
Table |I] is used to convert to log-plot, (b) The correspond- 
ing modified Avrami plot. The straight lines are obtained by 
least-square fitting. 



e. In contrast to the previous report for the Ising sys- 
tem [13] , a fortunate cancellation of the first and the sec- 
ond term of Eq. ([7]) does not occur and we have rather 
long incubation times in our phase- field model (Table 

Usually the exponential Arrhenius-type temperature 
dependence [27l. [28j 



tine oc exp 



(30) 



is assumed to analyze experimental incubation time. 
However, the incubation times obtained in TableUdo not 
increase exponentially but they increase linearly in Fig. [5] 
. The curve shows a nearly linear behavior expected from 
Eq.dUl): 



tu 



1 

— ; C>C — 

mk^T T 



(31) 



as the function of the inverse temperature 1/t cx 1/T 
with a slope that is inversely related to the undercooling 



e as the energy barrier W, is inversely proportional to e 
from Eq. ([221) [16]. Therefore when the undercooling is 
low (e = 0.2), the slope is steep and the incubation time 
increases rapidly as the temperature is lowered. Our nu- 
merical results in Fig. [5] seems to sug gest that the theory 
of Shncidman and Weinberg theory [la ] is qualitatively 
correct in our phase field model. 
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FIG. 5; The incubation time tine as the function of the 
inverse temperature 1/t ot 1/T. The incubation times do 
not increase exponentially but they increase linearly, and are 
roughly proportional to 1/t as predicted from Shneidman- 
Weinberg theory [l6l]in Eq. (|13p . The straight lines are ob- 
tained by least-square fitting. 

Figure [S| shows the temperature dependence of the 
growth time ig^. As intuitively expected, the growth time 
is an increasing function of the inverse temperature r^^. 
It increases as the temperature is lowered because nu- 
cleation rate becomes low. We should note that in our 
phase-field model, the interfacial velocity vq ~ v does not 
depend on the temperature from Eq. (I23p . Therefore, the 
temperature dependence of tgr in Eq. pop comes mainly 
from the temperature dependence of the nuclcation rate 

In fact, from Eq. pO[) we expect an activation- type 
temperature dependence 



(igr) OC 7^ (X exp 



knT 



(32) 



if the classical nucleation theory (CNT) is valid. There- 
fore, the Arrhenius plot of the growth time (Tgr) ver- 
sus inverse temperature should follow a straight 
line. Figure [7[ shows the Arrhenius (semi- log) plot of the 
growth time (Tgr)™. Both Tgr and m are those obtained 
from the least-square fitting of Eq. ([8[) and are tabulated 
in the first two columns of Table [H The curves are almost 
straight lines expected from CNT in Eq. ([32[l for high un- 
dercoolings e = 0.25 and e = 0.30. The classical picture 
of activation-type nucleation rate of CNT seems valid for 
high undercooling in our phase-field model. Then the to- 
tal time scale of the phase transformation ^inc+^gr that is 
the sum of the incubation time tjnc and the growth time 
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FIG. 6: The growth time tgr as the function of the inverse 
temperature 1/r oc l/T. 



cutoff ipc = 0.2 to prevent the thermal fluctuation ^(t, n) 
and the extra nucleation in our simulation. This choice 
'0c — 0.2 may not be appropriate since the critical nu- 
cleus is large and the energy barrier is high for low un- 
dercooling e = 0.2. In fact, the radius of critical nucleus 
become = a j tJi ~ 0.8 from Eqs. and (HH) as 
cr = 1/24 ~ 0.04 12| and A/i = 0.2/4 = 0.05 when 
e = 0.2, and roug hly -nRl ~ 2 to 3 neighboring cells 
should be transformed simultaneously in order to form 
a critical nucleus that can continue to grow. Then the 
small change of cutoff ■i/'c greatly influence the nucleation 
process and may affect the phase transformation dynam- 
ics. Therefore the deviations of simulation data from 
the prediction of CNT in Fig. [7| when the undercooling 
e — 0.2 could be due to the inappropriate choice of cutoff 
parameter ■(/'c and/or the limitation of controlling nucle- 
ation event by a single cutoff parameter. 



tgr, will not be described by a simple formula. Neither 
Eq. ^ predicted from CNT nor Eq. ^ of Shneid- 
man and Weinberg theory [l^ alone cannot describe the 
temperature dependence of total transformation time. 
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FIG. 8: The Avrami exponent m as the function of the inverse 
temperature 1/r oc l/T. The Avrami exponent m is much 
larger when the cutoff i/Jc is not introduced (t/jc = oo) than 
when = 0.20. 



FIG. 7: The logarithm of the growth time In [(tgr)"] pro- 
portional to that of the nucleation rate In [/s] as the function 
of the inverse temperature 1/r oc l/T. The logarithm of the 
nucleation rate In [7s] are roughly proportional to 1/r in ac- 
cordance with the CNT. The straight lines are obtained by 
least-square fitting. 

The curve deviated, however, from a straight line for 
the low undercoohng e — 0.2 in Fig.[71 Therefore, our nu- 
merical results seem to indicate that the CNT in Eq. (fTT|) 
cannot be used for the phase-field model when the un- 
dercooling is low. Looking at the curve in Fig. [7] for 
e = 0.2 closely, we note that the curve start to devi- 
ate from a straight line when the inverse temperature 
becomes large > 30 when the Avrami exponent m 
also becomes larger than the ideal value 3 (Table |T| . In 
fact, our cell dynamics code can not simulate the nucle- 
ation process correctly for low undercooling. It should 
be noted that we have artificially introduced an ad hock 



When the cutoff parameter is not introduced (0c — 
cxd), the thermal fiuctuation £,{t,n) and the extra nucle- 
ation may occurs even in the region where the nucleation 
has already started. Then the nucleation rate should be 
enhanced. The nucleation rate increases as more and 
more materials are transformed. The nucleation rate be- 
comes the increasing function of time I{t) ~ /pi" with 
non-zero exponent a ^ 0, which will lead to a larger 
Avrami exponent m from Eq. ([9]). Actually the Avrami 
exponent m in TableUbecomes larger than the ideal value 
3 for low temperatures > 30 at the low undercool- 
ing e — 0.2. Therefore the artificial cutoff tpc alone may 
not control the nucleation event in our phase field model 
properly. 

Figure [5] compares the Avrami exponents m calculated 
when tpc = 0.2 with those calculated when ipc — oo (with- 
out cutoff). The Avrami exponent m becomes much 
larger than ideal value m — 3 when ■i/'c = oo as expected. 
However, since we can still deduce not only the Avrami 
exponent m but the incubation time tjnc and the growth 
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time igr, the KJMA picture of phase transformation rep- 
resented by Eq. ([5]) is still qualitatively valid even when 
^Jc — oo. Therefore, qualitative picture of nucleation and 
growth will be still valid even if we do not introduce cut- 
off Tpc- 

So far, such a thermal fluctuation is expected to initiate 
nucleation of patter formation in phase field model [13, 
[T3 | or cell dynamics [l9„ 2Q3 29] , and is uncritically used 
without any special care such as our cutoff in this paper. 
Our phase field simulation for the most basic process of 
nucleation and growth clearly showed that such a naive 
expectation is only qualitatively valid and may not lead 
to quantitatively correct description of the dynamics of 
phase transformation. 

V. CONCLUSION 

In this report, we used the cell dynamics method 21] to 
study the whole dynamics of phase transformation from 
the homogeneous nucleation to growth in a simple phase 
field model. We found that the Kolmogorov-Johnson- 
Mehl-Avrami (KJMA) scenario of nucleation and growth 
is correct in this phase-field model. Specifically, the evo- 
lution of the volume fraction of transformed volume fol- 
lows the KJMA formula. By fitting the KJMA formula 
to the simulation data, we could deduce not only the 
Avrami exponent m but the incubation time tjnc and the 
growth time igi-. So far as the present author knows, 
there has been virtually no detailed quantitative study 
of the incubation time. Therefore our report is probably 
the first quantitative report of the theoretical study of 
the incubation time tine- 

The Avrami exponent m is found to be very close 
to the ideal value m = 3. The incubation time 
tine is inversely proportional to the absolute tempera- 
ture and is qualitatively in accord with the theory of 
Shneidman- Weinberg [T^. This conclusion would be 
valid for real materials unless some activation-type tem- 
perature dependence comes in through, for example, 
diffusion- limited growth. In fact, many experimental 
workers [l^, [Hi successfully analyzed their experimen- 
tal data assuming exponential activation-type tempera- 
ture dependence of incubation time without paying much 
attention to the Shneidman- Weinberg theory. However, 
this could partly due to the ambiguity of the definition of 
experimentally determined incubation time that could be 
actually the growth time. In our study we have separated 
the incubation time and growth time unambiguously us- 
ing the KJMA formula for the volume fraction. It turns 
out that both the incubation time and the growth time 
are the same order of magnitude in the phase- field model. 
Therefore, similar clear discrimination between the incu- 



bation time and the growth time, if possible, would be 
necessary to analyze experimental data quantitatively. 

The temperature dependence of nucleation rate / de- 
duced from the growth time igr follows the classical nu- 
cleation theory except at low undercooling. The error 
at the low undercooling is caused from the difficulty of 
controlling nucleation event in phase field model. An ad 
hoc cutoff to prevent thermal fluctuation and extra 
nucleation in the region where the phase transformation 
has already started cannot work well for the low under- 
cooling. Therefore, a blind use of thermal fluctuation 
in phase field model and/or cell dynamics method needs 
caution. It may not give a quantitatively correct dynam- 
ics of phase transformation. Our numerical simulation 
has unambiguously shown that an inappropriate choice 
of this cutoff parameter will lead to quantitatively incor- 
rect Avrami exponent m. 

In conclusion, we have shown that the most basic 
dynamics of phase transformation in the phase field 
model follows the scenario of Kolmogorov-Johanson- 
Mehl-Avrami. Furthermore, we have demonstrated that 
the incubation time ijnc and the growth time show 
different temperature dependence. Therefore, the total 
phase transformation time tine + ^gr will not be described 
by a simple single formula. If, however, t^c is dominant, 
that is expected when the temperature is lowered, the 
temperature dependence of phase transformation time 
follows usual activation-type formula of CNT. On the 
other hand, if iinc is dominant, the phase transformation 
time is not activation-type but is inversely proportional 
to the absolute temperature [16]. 

Finall y, i n contrast to the standard cellular automa- 
ton [33, mi and the lattice model [H, |33|] where arti- 
ficial evolution processes are introduced algorithmically, 
the evolution in the phase-field model is driven by the free 
energy and the thermal noise, and is completely free from 
artificial parameters. Therefore, the phase-field model 
with cell dynamics method is powerful method to study 
the various phase transformation scenarios, in particular, 
of the time evolution including the incubation time and 
the growth time without introducing many unknown and 
uncontrollable parameters. For example, various scenar- 
ios of heterogeneous nucleation are already clarified [13] 
by this phase-field model with cell dynamics method. 
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